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We re-analyze the production of seed magnetic fields during Inflation in (R/m 2 )" F^F^ and 
IFfivF 11 " models, where n is a positive integer, R the Ricci scalar, m a mass parameter, and I <x r] a 
a power-law function of the conformal time r], with a a positive real number. If m is the electron 
mass, the produced fields are uninterestingly small for all n. Taking m as a free parameter we find 
that, for n > 2, the produced magnetic fields can be sufficiently strong in order to seed dynamo 
mechanism and then to explain galactic magnetism. For a > 2, there is always a window in the 
parameters defining Inflation such that the generated magnetic fields are astrophysically interesting. 
Moreover, if Inflation is (almost) de Sitter and the produced fields almost scale-invariant (a ~ 4), 
their intensity can be strong enough to directly explain the presence of microgauss galactic magnetic 
fields. 
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I. INTRODUCTION 

Microgauss magnetic fields are present in any type of 
galaxies Up to now, astronomical observations have 
not been able to clarify the nature and origin of these 
fields, so that "galactic magnetism" is still an open issue 
in cosmology. 

If on the one hand, the presence of magnetic fields in all 
galaxies could be explained by the action of astrophysical 
mechanisms operating during or after large-scale struc- 
ture formation 0], on the other hand, the detection of 
magnetic fields in galaxies at high redshifts Q represent 
a strong hint that magnetic fields could have been gener- 
ated in the very early Universe. Clearly, the detection of 
magnetic fields in the Cosmic Microwave Background ra- 
diation would prove the primordial origin of cosmic mag- 
netic fields Q. In such a case, it is very probable that 
they have been generated during an inflationary epoch 
of the Universe, since only in this case it is possible to 
generate fields correlated on comoving scale as large as 
galactic ones. 

Indeed, many mechanisms 0, @, 0, d, Q , acting during 
Inflation, have been proposed to produce seed fields since 
the seminal paper by Turner and Widrow [lj| • 

It is worth noting that, in order to explain galactic 
magnetism, one needs the presence of seed magnetic fields 
prior to galaxy formation of intensity generally much less 
than lfiG. In fact, when a protogalaxy collapses to form 
a galactic disk, magnetic fields suffer an amplification due 
to magnetic flux conservation. Moreover, due to magne- 
tohydrodynamic turbulence effects and differential rota- 
tion of galaxy, seed fields can be further amplified. This 
last mechanism, know as "galactic dynamo" ll|, can 



be very efficient and, in principle, allow extremely week 
seeds to reproduce the properties of presently-observed 
galactic fields. 

If dynamo is operating, the present-day intensity and 



correlation length required for a successful amplification 
are 1 B > 10~ 33 G and A > lOkpc while, without dy- 
namo, a comoving seed field B > 10~ 14 G correlated on 
comoving scales of order IMpc is needed to explain galac- 
tic magnetism Note that, following Ref. we are 
assuming that the redshift at which a protogalaxy sepa- 
rates from the Hubble flow to then collapse is z pg ~ 50. 

In this paper, we re-analyze the generation of seed 
magnetic fields during Inflation considering non-standard 
terms in the electromagnetic Lagrangian of the form 
(R/m 2 ) n F^F^ and IF ftu F>" / , where n is a positive 
integer, R the Ricci scalar, m a mass parameter, and 
I{4>, R, ...) a scalar function depending on any scalar field, 
gravity, and on any other background field. 

The first model was first proposed by Turner and 
Widrow for the case n = 1 and then extended to the 
case n > 1 by Mazzitelli and Spedalieri Q. They found 
that, taking m = m e (m e being the electron mass), astro- 
physically interesting seed fields can be generated during 
Inflation for adequate values of n. In Section II, we will 
find a different result from Ref. jfj. Indeed, if m = m e , 
the produced fields are uninterestingly small for all n. 
Nevertheless, taking m as a free parameter, the possibil- 
ity to generate magnetic fields sufficiently strong in order 
to seed dynamo mechanism is not excluded. 

The second model was first investigated by Bamba and 
Sasaki [3] who assumed de Sitter Inflation and I to be 
a power-law function of the conformal time T), I oc r) a , 
with a a negative real number. In Section III, we will 
extend this model to the case of positive real a for both 
Power-Law and de Sitter Inflation. 



1 It is worth noting that galactic dynamo is a still-debated and 
model-dependent amplification mechanism, and that the mini- 



mum value of the seed field used in this paper, B 



10" 



G, 



corresponds indeed to an extremely efficient dynamo action. For 
details see, e.g., Ref. 0] and references therein. 



2 



In Section IV we draw our conclusions. 



II. R n F 2 MODELS 



We consider a photon described by the effective 
Lagrangian 



r — __ F F^ u — - TF F^ v 
I- — ^ r ^ v r 1 r^ u r 



—2) ID 



where F„ v — d^A v — d v A„ is the electromagnetic field 
strength tensor, R the Ricci scalar, m a parameter having 
dimension of a mass, and n a positive integer. (From now 
on, Greek indices range from to 3, while Latin ones from 
1 to 3.) 

We will work in a flat Universe described by a 
Robertson- Walker metric ds 2 — a 2 (dr] 2 — rfx 2 ), where the 
conformal time r\ is related to the cosmic time t through 
dr\ = dt/a. We normalize the expansion parameter so 
that at the present time t , a(t ) = 1. 

We consider the case of Power-Law Inflation de- 
scribed by the equation of state ptot = 7Ptot with 
— 1 < 7 < — 1/3. In this case, the Hubble param- 
eter H 2 = (87r/3)yO t ot/m-p 1 and the Ricci scalar R = 
87r(l — 37)ptot/m-pi evolve as H 2 cx R cx ptot °c if, 
where s = -6(1 + 7)/(l + 37). Here, p tot (p tot ) is 
the total energy (pressure) density during Inflation and 
mpi ~ 1.22 x 10 19 GeV is the Planck mass. 

We assume that during Inflation R ^> m 2 , so that we 
can neglect the Maxwell term in Lagrangian ([1]). Since 
R is a decreasing function of time, it suffices to impose 
that i?i < m 2 or equivalently m < m max , with 



VM1-37) (—), (2) 



mpi 



where M 4 = ptot(?7i) is the total energy density at the 
end of Inflation. Here and in the following, the subscript 
"1" indicates the time when Inflation ends. 

The spectrum of gravitational waves generated at In- 
flation is submitted to constraints coming from Cosmic 
Microwave Background analysis which requires p tot to be 



less than about 10 



on the scale of the present Hub- 



ble radius [Io| . This turns in a upper limit on the value 

< s max , where [13] 
f log 10 (M/mpi)]. 



of M and s: 



M < lCT 2 m P1 and s 
41og 10 (M/m P1 )]/[29.8 



One must also impose that M > lGeV, so that the pre- 
dictions of Big Bang Nucleosynthesis are not spoiled [l0( . 

We work in the Coulomb gauge, A = J^j diA{ = 0, 
and we expand the electromagnetic field — (Aq , A) 



A(t7,x) 



d 3 k 



(27r)3V2fc a tr 



where k = |k| and £k,a are the transverse po- 
larization vectors satisfying the completeness relation 

Ea( £ k,aM£k,Ji = ^y( k )> wi t h PijQ*) = &ij - kikj/k 2 . 
In order that the annihilation and creation operators 
Ok,a and a k satisfy the usual commutation relations 
[ak lCt , a '] — (27r) 3 (5 QQ '(5(k — k'), we must impose the 
normalization condition Ak (f?)-^ (rf) — A* k (rf)Ak(rj) = i/I, 
where a dot denotes differentiation with respect to rj. 
This assures, in turn, that the field A^ and its canoni- 
cal conjugate momenta ttq = and iti = IAi satisfy the 
usual commutation relations: 



[Ai(77,x), 7^(77, y)] = i 



d 3 k 
(2tt)3 



s <k(x-y)p ii(k)i (4) 



The equation of motion for A k {rf) follows from La- 
grangian (J): 

A k ( V ) + — A k ( V )+k 2 A k (r 1 )=0. (5) 
77 

The solution of the above equation, taking into account 
the normalization condition on A k (rf), is easily found: 



r(l+2,)/4 7 -l/2 1^11/2^(2) (g) 



where = (ns — l)/2 and if,) (x) is the Hankel function 
of second kind. (Observe that, since both n and s are 
positive, it results v > —1/2). 

We are interested in the study of large-scale elec- 
tromagnetic fields, that is in modes whose physical 
wavelength is much greater than the Hubble radius, 
Aphys 3> H^ 1 or equivalently \krj\ <C 2/|l + 37 1, where 
Aphys = o,X, A = 1/k is the comoving wavelength, and 
we used ry = 2/[(l + 3j)aH], Therefore, in Eq. ©, 
we can replace the Hankel function with its small- 
argument expansion: H^\x) ~ i-K~ lr L v Y{y)x~ v for 
1/ > 0, h¥\x) ~ 2i7r- 1 lna; for 2/ = 0, and i^ 2) (ir) ~ 
-m- x 2- v Y{-v)e-" v x v for < 0, where T(x) is the 
Gamma function. 

After Inflation, the Universe enters in the so-called re- 
heating phase, during which the energy of the inflaton 
is converted into ordinary matter. In this paper, we re- 
strict our analysis to the case of instantaneous reheating, 
that is after Inflation the Universe enters the radiation 
dominated era. In radiation era, the general expression 
for the electromagnetic field is 



ikrj 



(7) 



(3) 



where at and (3 k are the so-called Bogoliubov coeffi- 
cients [12j , determining the spectral number distribution 
of particles produced from the vacuum. By matching 
expressions © and Q and their first derivative at the 
time of the end of Inflation, we find the spectrum of the 
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electromagnetic field generated from the vacuum at large 
scales: 



\Ar(vi)\ 2 = \Pk\ 



[2T(1 + v)} 2 
2tt7i \kr ll \ 1 + 2v 



(8) 



valid for 1/ > 0. When — 1/2 < f < 0, the electromagnetic 
vacuum fluctuations go like |A£ ac (?7i)| 2 <x |fc?7i| 1+2!y and 
then are vanishingly small for \kr}x\ C 1. In the case 
v = —1/2, conformal invariance is recovered and then, 
as it should be, we find |^4fc ac (f?i)| 2 = 0, exactly, for all 
\k m \. ' ' 

We can define, now, the average magnetic field on a 
comoving scale A as @ 



B 2 x (v)- 



d 3 yW x (\x-y\)VxA(r,,y)\ 2 \0), (9) 



where Wa(|x|) is a suitable (real) window function, and 
the vacuum state |0), such that ak, Q |0) = 0, is normalized 
as (0 1 0) = 1. Taking into account Eqs. ([3]) and ([9]), we 
obtain 



(10) 



where W\(k) is the Fourier Transform of the window 
function and 



11) 



2n 2 a 4 

is the magnetic power spectrum. As windows func- 
tion we can take a gaussian window, IFa(|x|) = 
(2^A 2 )- 3 / 2 e-l x l 2 /( 2A2 \ so that W x {k) = ' 2 . 

From the end of Inflation until today, due to the 
high conductivity of the cosmic plasma, the magnetic 
field evolves adiabatically [10|, a 2 B\ = const, so that 
^today = a 2 Bx ( m y Inserting Eq. p in Eq. $11$ and 
taking into account Eq. (fT0|) . we finally obtain: 



B 



h 



today 



2 , T(l + ^)[r(3/2- J /)] 1 /2 



(2tt) 3 / 2 I, 



1/2 



-1/2 



m J 



V3 



27r|l + 3 7 | 



J*s 



m) 



g*s(To) 



1/3 



MTo 



(12) 
(13) 
(14) 



In Eq. (|T4]) , we used the fact that during radiation and 
matter dominated eras the expansion parameter evolves 
as a cx g^g^T^ 1 , where T is the temperature and g*s{T) 
the number of effectively massless degrees of freedom re- 
ferring to the entropy density of the Universe [HI ■ The 
temperature at the end of Inflation, the so-called reheat- 
ing temperature, is [ToJ T\ = M, while To is the actual 
temperature. 2 
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FIG. 1: Actual magnetic field as a function of M for m = m e 
(upper panel) and m = m max (lower panel) at the comoving 
scale A = lOkpc for different values of n: n = 1 (continuous 
line), n = 2 (long-dashed line), n = 3 (dashed line), n = 4 
(dot-dashed line), n = 5 (dotted line). The horizontal dotted 
lines refers to the minimum seed field required for dynamo 
amplification, B ~ 10 _33 G. In the lower panel, on the left 
(right) side of the vertical dotted line it results m < m e (m > 
m e ). 



It is important to observe that, in order to avoid an 
infrared divergence in the magnetic power spectrum, we 
must impose that v < 3/2. 

The "best case scenario", i.e. the case of maximum 
strength of the magnetic field, corresponds to s = s max . 

In Fig. 1, we show the actual magnetic field, in the best 
case scenario, as a function of the Inflation energy scale 
M for different values of the index n at the comoving 
scale A = lOkpc in the two cases m — m e (upper panel) 
and m = m max (lower panel). The range of M is such 
that < v < 1.49 (both panels) and m < m max (upper 
panel) . 

As it is clear from the figure, if m = m e the produced 
fields are not astrophysically interesting. This disagrees 
with the results discussed by Mazzitelli and Spedalieri in 



In the following we will use the values [IS 



T 



2.35 x 



10- 13 GeV, g tS (To) ^ 3.91, and g* s (Ti) = 106.75 (referring 
to the number of effectively massless degrees of freedom of Stan- 
dard Model). It is useful to recall that 1G ~ 6.9 X lCT 20 GeV 2 
and IMpc ~ 1.56 x 10 38 GeV _1 . 
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Ref. @ where it is argued that, for adequate values of n 
and 7, the produced fields can serve as seed for the cos- 
mological magnetic field. The discrepancy resides in the 
fact that, the astrophysically interesting fields they ob- 
tained correspond, indeed, to values of n and 7 (and then 
of v) such that the related power spectrum is infrared- 
divergent, and then are non-physical. 

In the "optimistic" case m = m max , for all n > 2 
there is a range in the parameter M such that the inten- 
sity of the produced magnetic field is greater than the 
minimum seed field required for dynamo amplification, 
B > 1(T 33 G. 

III. IF 2 MODELS 

We consider now a model given by Lagrangian (fTJ) with 
/ an arbitrary power function of the conformal time: 




where I\ is a constant and we assume that a is a positive 
real number. The above parametrization of the function 
/, already used in the literature Q, is not the unique 
admissible one. However, leaving I an arbitrary function 
of the conformal time, if on the one hand would render 
our analysis more general, on the other hand would allow 
us just to get qualitative results. 

We study the production of seed fields during an in- 
flationary epoch described by both Power-Law and dc 
Sitter Inflation. In the latter case, the equation of state 
describing the evolution of the Universe is ptot = 7Ptot 
with 7 = — 1, so that the Hubble parameter is a constant 
and 77 = — l/(aH). 

If 1% > 1, we can neglect the Maxwell term in Eq. (fT]) 
and then the analysis performed in Section II applies also 
to the case at hand. In particular, Eqs. (0, ©, (O, an d 
(I14p are still valid provided ns is replaced by a. Conse- 
quently, now it results v = (a — l)/2 and the condition 
< v < 3/2 translates to 1 < a < 4. The expression for 
the actual magnetic field can be recast as: 

R toda y _ 10 12 - 55 [r(2-a/2)]V 2 /Jtf V*' 2 (q _ 4)/2 
CaJ i/2 |l + 3 7 |-«/ 2 UpJ 10kpC 

(16) 

where Aiokpc = A/(10kpc) and ( a is an increasing func- 
tion of a of order unity (such that £i ~ 0.6 and (4 ~ 2.3). 

In Fig. 2, we show the value of actual magnetic field 
both for Power-Law (in the best case scenario corre- 
sponding to s = s max ) and de Sitter Inflation as a func- 
tion of M (which in this case can assume all values in 
the interval lGeV < M < 10 _2 mpi) for different values 
of the parameter a at the comoving scale A = lOkpc and 
for h = 1. 

We find that for a > 2 there is always a minimum 
value of M, M m i n , such that if M > M m i n the generated 
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FIG. 2: Actual magnetic field in the case of Power-Law (con- 
tinuous lines) and de Sitter (dashed lines) Inflation as a func- 
tion of M for different values of a. From down to top: 
a = 1, 2, 3, 3.99. The first two curves from the top correspond 
to A = IMpc, while the other ones to A = lOkpc. The horizon- 
tal dotted lines refer to the minimum seed fields required for 
dynamo amplification, B ~ 10 _33 G, and to directly explain 
galactic magnetism, B ~ 10~ 14 G. 



magnetic field is greater than the minimum seed field re- 
quired for a successful dynamo amplification. Moreover, 
if the field is almost scale- invariant, i.e. a ~ 4, and if 
M is large enough, its intensity can be strong enough 
to directly explain galactic magnetism (see the first two 
curves from the top in Fig. 2 referring to A = IMpc). In 
fact, expanding Eq. (fl"6|) about a = 4 and taking 7 ~ — 1 
(which corresponds to take large values of M), we find 

B?** 2Xl(r7 f^TG, (17) 

giving Bf day > 1(T 14 G for M > 3(1 - a/4) 1 / 4 10 15 GeV. 

III. CONCLUSIONS 

Why are large-scale magnetic fields present in all galax- 
ies is still a debated question in cosmology. They ori- 
gin could reside in an amplification of Inflation-produced 
electromagnetic vacuum fluctuations through galactic- 
dynamo mechanisms. 

In this paper, we have re-analyzed the production of 
seed magnetic fields during Inflation considering not- 
conformal-invariant terms in the electromagnetic La- 
grangian of the form {R/m 2 ) n F^F^ and IF^F^. 

The first model was first proposed by Turner and 
Widrow [l(| for the case n = 1 and then extended to 
the case n > 1 by Mazzitelli and Spedalieri [6j. 

The second model was first investigated by Bamba and 
Sasaki [7] in the framework of de Sitter Inflation and 
assuming / to be a power-law function of the conformal 
time 77, / oc 77", with a a negative real number. 

Concerning the first model, Mazzitelli and 
Spedalieri [6j found that, taking m — m e (m e be- 
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ing the electron mass), astrophysically interesting seed 
fields can be generated during Inflation for adequate 
values of n. Our analysis instead, as discussed in Section 
II, shows that if m is the electron mass, the produced 
fields are uninterestingly small for all n. Nevertheless, 
taking m as a free parameter, the possibility to generate 
magnetic fields sufficiently strong in order to seed 
dynamo mechanism is not excluded. 

Regarding the second model, we have extended it to 
the case of positive real a for both Power-Law and de 
Sitter Inflation. We have found that, for a > 2, there is 
always a window in the parameters defining Inflation such 
that the generated magnetic fields are astrophysically in- 
teresting. Moreover, if Inflation is (almost) de Sitter and 
the produced fields almost scale-invariant (a ~ 4), their 
intensity can be strong enough to directly explain the 
presence of microgauss galactic magnetic fields. 
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